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Abstract- A generalized algorithm for multiplication is 
proposed through recursive application of the Nikhilam Sutra 
from Vedic Mathematics, operating in radix - 2 number 
system environment suitable for digital platforms. Statistical 
analysis has been carried out based on the number of 
recursions profile as a function of the smaller multiplicand. 
The proposed algorithm is efficient for smaller multiplicands 
as well, unlike most of the asymptotically fast algorithms. 
Further, a basic block schematic of Hardware Implementation 
of our algorithm is suggested to exploit parallelism and speed 
up the implementation of the algorithm in a multiprocessor 
environment. 
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I. Introduction 

Vedic Matliematics: 

• Vedic Mathematics deals mainly with 16 Sūtras and 
their applications for carrying out tedious and 
cumbersome arithmetical operations, and to a very 
large extent, executing them mentally; see[5]. 

• Nikhilam Navatascaramam Dasatali Sūtra (or simply 
Nikkilam Sūtra) is one of these 16 Sūtras used for 
multiplication and has been successfully applied to 
overcome drawbacks of conventional schemes. 

This is where the fields of modern computing and Vedic 
Mathematics converge. 

Need For Efficient Multiplication Algoritlmi: 

• Use of numerical methods is prevalent in most 
software algorithms. Such applications demand an 
efficient code for basic mathematical operations, one 
of them being multiplication. 

• Real Time Systems demand instantaneous response to 
environmental variables and quick execution of taken 
decision. 

• Multiplication algorithms find applications in Digital 
Signal Processing fDSP) for discrete Fourier 
transforms, Fast Fourier transforms, convolution, 
digital filters, etc. Therefore any new multiplication 
algorithm opens up a new approaclt for improving 
existing schemes. 


Existing Algorithms: 

• The naive binary multiplication algorithm has a time 
complexity Ofn A 2) where n is the number of bits of 
the numbers being multiplied. 

• The Schōnhage-Strassen algoritfun is an 

asymptotically fast multiplication algorithm for large 
integers. It was developed by Arnold Schōnhage and 
Volker Strassen in 1971. The run-time bit complexity 
is, in Big O notation, Ofn log n flog log n)), while the 
arithmetic complexity is Ofn log n). The algorithm 
uses recursive Fast Fourier transforms in rings with 
elements, a specific type of number theoretic 
transform. 

• The Schōnhage-Strassen algorithm was the 

asymptotically fastest multiplication method known 
from 1971 to 2007 when a new method, Furer's 
algorithm, was announced with lower asymptotic 
complexity. 

• The drawback of the Schōnhage-Strassen algorithm is 
that it starts to outperform older methods such as 
Karatsuba and Toom-Cook multiplication for 
numbers beyond 10,000 to 40,000 decimal digits and 
so is the case with most of the asymptotically fast 
algorithms. Hence these algorithms are not suitable 
for general purpose digital systems. 

• There is no algoritfun in our knowledge whose 
complexity depends upon ratio of Os and 1 s in 
multiplicands rather than magnitude of multiplicand. 

This motivated us to come up with an algoritfun which is 
‘time efficient’ over the entire range of numbers for more 
generalized applications. 


II. Explanation Of Nikhilam Sūtra 

One of the 16 Sūtras of Vedic Mathematics, Nikhilam Sūtra 
stated algebraically as follows: 

Consider two numbers and n 2 such that 

n r — (x — a) a — (x — n r ) 

n 2 — (x — b) b — (x — n 2 ) 

where x = base, a,b — differences from base, then 

n t x n 2 = ( x - a )( x - b ) = x((x - a) + (x - b)- x ) + a* b 


This calls for a ‘time efficient’ algorithm for Thus the Nikhilam Sutra effectively breaks up a large 

‘multiplication’ to improve processor throughput. multiplication n t X n 2 into a small multiplication (a X b) 

and addition | (x — a) + (x — b) — x]. 



Case 1.1: Base as power of 10, «, and n 2 less than base: 
Consider n t = 99 and n 2 = 98 


Case 2: Modified Base, Multiplication Method: 
Again consider n t =49 and n 2 = 48 


EXAMPLE 1.1 


Number 

Difference 

Base x 

O-a) 

a 


99 

+01 


(x-b) 

b 

100 

98 

+02 


(x - a) + (x - b)- x 

97 

ab 

02 



Hence answer = 9702 


EXAMPLE 2 


Number 

Difference 

Base x 

(x-a) 

a 


49 

+i 


(x-b) 

b 

10*5= 50 

48 

+2 


{(x - a) + (x - b)-x) * 5 

235 

ab 

2 



Hence answer = 2352 


Case 1.2: Base as power of 10, tij and n 2 greater than base: 
Considern-L = 102 and n 2 = 101 


EXAMPLE 1.2 


Number 

Difference 

Base x 

( x — a) 

a 


102 

-02 


(x-b) 

b 

100 

101 

01 


[x - a) + (x - b) - x 

103 

ab 

02 



Hence answer = 10302 


Case 1.3: Base as power of 10, base between n 2 and n 2 : 
Consider n, = 101 and n 2 = 99 


EXAMPLE 1.3 


Number 

Difference 

Base x 

(x-a) 

a 


101 

-01 


(x-b) 

b 

100 

99 

+01 


(x - a) + (x - b) - x 

100 

ab 

-01 



Note: 

• The Modified Base Method is applied for all three 
conditions of n 2 and n 2 as in Cases 1.1, 1.2 and 1.3 
using either Multiplication or Division shown above. 

• If one of the numbers to be multiplied is a Base itself, 
then we have to consider only the calculation on LHS, 
as product of differences is zero. 

• Number of digits for product (a X b) equals the power 
of 10 with which base is associated. 


Extension to Radix-2 Number System: 

The Nikhilam Sūtra as described in [5] is applied 
to decimal, i.e., Radix - 10 system only. We have extended 
the domain of Nikkilam Sūtra to encompass Binary, i.e., 
Radix - 2 number system. This enables us to harness it’s 
prowess for implementing the unsigned Multiplication 
algorithm on digital platform. 

All cases discussed for decimal system are applied 
to hexadecimal in analogous manner and are found to work 
accurately. 

The implementation can also be easily extended to 
perform signed binary number multiplication. This can be 
accomplished by processing the MSBs i.e. sign bits of both 
multiplicands separately. If ‘1’ denotes positive and ‘0’ 
denotes negative number, then the sign of the 
Multiplication result can be obtained by simple EX-OR 
logic applied on the sign bits of the two numbers. 


Sign Bitofnl 
Sign Bit of n2 


Sign Bit of (nl x n2) 



The product ab being negative is to be subtracted from 
10000. Hence answer = 9999 


























III. Need For Recursion 


IV. Recursive Algorithm For Binary Multiplication 


Review ofExisting Works: 

The comparison with the existing similar attempts 
mentioned in references leads to the following 
observations: 

• Our previous work “An Efficient Binary 
Multiplication Algorithm based on Vedic 
Mathematics” see [1] is an 8 bit non recursive binary 
multiplication method. The 8 bit binary multiplication 
method in [1] is much superior to that in [3]. Even in 
the worst case the former is faster by over 50% and 
much faster for the average case. 

• The method in [2] compares the Nikbilam Sutra and 
Urdha Tiryak Sutra and highlights the cases when 
either of the algorithms is more efficient than the 
other. Further the paper explains the basic mechanism 
of Nikhilam Sutra applied to digital arithmetic but is 
restricted to 4 bit binary multiplication. For higher bit 
multiplication the algorithm is not extendible easily 
and needs to consider a very very large number of 
cases of multiplicands and hence it is unacceptable 
when compared to [1]. 

• The ‘efficient’ multiplication algorithms like 
Schonhage-Strassen algorithm, Fiirer's algorithm only 
perform well for very large numbers as multiplicands 
and hence not suitable for general purpose 
applications. 

• The Algorithm in [1] can be extended to higher bit 
multiplications by breaking up the number into 8 bit 
segments and operating separately on them however 
the proposed generalization method is much more 
efficient than such an implementation. 


Scopefor Recursion: 

The Nikhilam Sutra based multiplication can be 
summarized as the expression of a large multiplication into 
a small multiplication and an addition. For the Modified 
Base method, this involves another small multiplication for 
base adjustment. 

The proposed method completely avoids using the 
Modified Base method by appropriate base selection and 
hence avoids additional small multiplication. The possible 
bases in the proposed method are in geometric progression 
as opposed to arithmetic progression in [1]. 

Due to such base selection, the differences from 
Base may not be small as expected by Nikhilam Sutra to be 
efficient. Thereby in such case, the product of differences 
may itself be performed by recursively applying the 
proposed algorithm. This process continues until the final 
product of differences is very small so as to be 
accomplished with multiple additions. 

In this manner any two numbers, irrespective of 
their magnitude, can be decomposed successively by the 
proposed algorithm. Effectively the multiplication is broken 
down into series of additions, as explained in the following 
section with examples. 


We propose a novel algorithm for binary 
multiplication which depends on the ratio of the number of 
ls and the number of Os in the binary representation of the 
highest of the two multiplicands. 

In this section we first introduce the algorithm, and 
then explain it with an example. The next section provides a 
theoretical proof based on the principle of Mathematical 
Induction. 

Algoritlim: 

1. Of the two multiplicands, say X and Y, select the 
smaller multiplicand. 

2. If either X or Y is 0 or 1 then go to step 6. 

3. If the second MSB of this smaller multiplicand is 1 then 
select ‘Base’ B as 2 A b else the select B as 2 A (b-l). 

4. Subtract this base from both the multiplicands to get 
differences X' and Y'. 

5. Apply steps 1 to 3 to X’ and Y’ recursively till X’ = 1 or 
Y' = 1 OR X' = 0 or Y' = 0. 

6. At recursions, get the partial result at each step by either 
X + Y' OR Y + X'. 

7. Multiply X’ and Y’ in the final stage (where one of X' or 
Y' is either 0 or 1. 

8. At each step the respective b no. of bits of the partial 
result are considered while other MSBs are carries 
generated. 

9. Using step 8, merge the partial results to get the final 
product. 


Example: 

Consider the following example illustrated in 
binary and verified in the decimal number system. Each 
column stands for a recursion and shows the partial results 
obtained at that recursion. The Multiplication result, after 
merging the partial results is shown in penultimate row, and 
verification in decimal system is shown in the last row. 

23 X 21 


Base = 16 
10111 - 10000=111 


Base = 4 



10111 

111 

11 

10101 

101 

01 

101100 

11100 

1000 

11 


00 

11 

10 



11110 



111100011 

23 X 21 = 483 
















V. Proof By Mathematical Induction 


Step 3 : Proof for base 2 6,1+1 using assumption in 
previous case 



Therefore, Assumption becomes : 


+ 


rij + 

(n 2 - 2 61 ) 

2 &! 

2 b i - 2 fc 2 



(% - 2 bl ) + 

(n 2 - 2 &1 — 2 &2 ) 

2*2 

,..-2 6 fc) x 


+ 


( n l ■ 

(n 2 - 2 61 - 2 6 2 - .- 2 6fc ) 


= n, x n 2 


ī.e. 

% + 

(n 2 — 2 61 ) 
2 61 



(% - 2 61 ) + 

( n 2 — 2 6l -2 6z ) 


+ 


= (% x n 2 ) - 

(ni - 2 61 - 2 &2 - .- 2 6fc ) x 

(n 2 - 2 61 - 2 62 - .- 2 6fc ) 


.Equation (1) 


Using result from equation (1) we get: 

, . Yni - 2 61 - 2 62 - .- 2 6fc ) x 

= (n, x n 2 ) - . . . 

(n 2 - 2 bl - 2 62 - .- 2 6fc ) 

4cA-.,H 

( 2 6fc+1 ) 

+ [(x- 2 6fc+1 ) x ( y — 2 6fc+1 )] 

= (n x x n 2 ) — [x xy] + |[ (y — 2 6fc+1 )] X | 

( 2 6fc+1 ) 

+ [(x-2 6fc+1 ) x (y — 2 6fc+1 )] 

= (% x n 2 ) - [x x y ] 

+ x 2 6fc+1 + y x 2 6fc+1 
_ 2 26fc+1 } 

+ [x xy- x x 2 6fc+1 - y x 2 6fc+1 
+ 2 26fc+1 ] 

= niXn 2 -xxy + xx 2 6fc+1 + y x 2 6fc+1 - 2 26fc+1 
+ x x y - x x 2 6fc+1 - y x 2 6fc+1 

-I- 2 26fc+1 

= n. x n. 


Thus above Induction proves the correctness of the 
algorithm and it holds for all bases. 






















VI. Hardware Implementation 


VII. Statistical Analysis 


Determining the base: 



<D 

CZ3 

03 

PQ 


T3 

8 


3 

<u 


04 


The base calculator allows determination of required base 
in one clock cycle. 

• Consider the input operand 00110110 

• The priority encoder output will be A 2 A X A 0 — 101 

• Effectively the decoder will output the required 
base viz. 00100000 


Multiplier Hardware: 






The plots describe the number of recursions 
required for each multiplicand from 0 to 127, where that 
multiplicand is the smaller of the two multiplicands. X-axis 
represents the smaller multiplicand ‘n’ and Y-axis shows 
corresponding number of recursions. Intervals shown on x- 
axis are the points where the number of bits ‘b’ changes, 
where ‘b’ is the number of bits needed to represent ‘n’. 

It can be easily seen form the plots that there is a 
pattern in the recursions for every ‘b’ bit number. The 
information about the number of occurrences of worst case 
i.e. number of maximum recursions cases for any 
multiplication involving smaller multiplicand of ‘b’ bits can 
be seen from the plots. 

If f(b) represents this count of cases with 
maximum number of recursions required for any 
multiplication involving smaller multiplicand of ‘b’ bits, 
then we have the recurrence relation for f(b) from the plots 
as follows. 

Initial Values: 

f(l) = 0, f(2) = 0, f(3) = 1, f(4) = 5 
Recurrence Relations: 



= base calculator 


If ‘b’ is odd and b>3 

f(b) = 2*f(b-2) 

If ‘b’ is even and b>4 

f(b) = 2*f(b-2) + f(b-l) 


where, 
























































































































































































From the above relations, f(b) can be calculated 
for any ‘b’ which gives us quantitative idea about number 
of worst case occurrences for each ‘b’. 

The equations show that for odd values of ‘b’, the 
cardinality of worst case recursions increases by a fixed 
multiple of 2 over that of ‘b-2’ where as the corpus and the 
range increase as a power of 2. For the case where 'b' is 
even valued, the cardinality of worst case recursions 
increases only linearly depending on cardinalities for ‘b-l’ 
and ‘b-2’, whereas the range still increases in powers of 2. 

The worst case recursions being ceil(b/2) the worst 
case complexity of the algorithm will be 0(n*b). 

VII. CONCLUSION AND FUTURE WORK 

We presented an algorithm which can be used in 
future to compute multiplication of two variable bit 
numbers. The algorithm was found to solely depend on the 
ratio of the number of l's and O's used to represent a 
number in binary, rather than on the magnitude of the 
number. It was seen frorn the graphs that as the ratio 
approaches 1 the number of operations required for the 
multiplication increases and decreases as the ratio tends to 
move close to 0. The statistical analysis provided a clear 
idea of the efficiency of the algorithm over the range of 
input multiplicands. The algorithm was found to be 
implementable by the hardware proposed and is much more 
suitable for applications on multiprocessor environments, 
than those proposed in [1], [2] and [3]. 

In future we plan to work on the average case 
analysis of the algorithm we propose. We intend to do the 
hardware simulation of the algorithm to analyze the 
empirical results of the multiplications over a corpus of 
numbers. 
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